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ABSTRACT 

X-ray spectroscopy of cooling-flow clusters reveal an unexpected deficit of soft X-ray 
emission lines from gas at ~ 1/3 of the ambient plasma temperature, across a wide 
range of X-ray luminosities and virial temperatures. We propose excess ionization from 
either a population of suprathermal electrons or photoionization by X-ray continuum 
emission from hot gas or a central AGN as a means of decoupling the thermal state 
of the gas from its emission line spectrum. The former effect is thought to operate 
in the solar corona. Because they generally become important at some fixed fraction 
of the cluster gas temperature, such mechanisms could in principle explain both the 
universality and temperature dependence of the emission line suppression, properties 
which none of the present-day models based on gas heating can explain. Ultimately 
these models cannot explain the observations; however, they have attractive and robust 
features which could be useful in elucidating a final solution to the soft X-ray deficit. 



1 INTRODUCTION 



Cooling flow clusters have long been a mystery. While the 
X-ray luminosity is consistent with gas cooling rates of up 
to 100 — IOOOM0 yr~ , whether such a cooling flow is indeed 
really operating and the fate of the cooled gas has been a 
matter of long-standing controversy (see Fabian (1994) for a 
review) . X-ray spectroscopy acquired by the Reflection Grat- 
ing Spectrometer on XMM-Newton has added yet another 
puzzle to the picture: there is a strong deficit of emission 
lines associated with the low temperature gas which should 
be present in a cooling flow. In particular, the observed dif- 
ferential emission measure from clusters is consistent with 
(Peterson et al 2002): 
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where M is the inferred mass dropout rate, T is the back- 
ground temperature, and a ~ 1 — 2. This conflicts with the 
standard isobaric cooling model, which predicts a = 0, inde- 
pendent of the details of gas clumping or geometry (Fabian 
(1994) and references therein). The observed emission mea- 
sure of low temperature gas is therefore lower than expected. 
In particular, while the spectra exhibits strong emission from 
the ambient temperature T a down to about T D /2-which sug- 
gests that cooling is indeed taking place-there is signifiantly 
less emission from lower temperature gas than is expected. 
What is particularly striking is that this seems to be a uni- 
versal phenomenon, across a wide range of cluster tempera- 
tures and mass deposition rates: emission from gas at ~ T /3 
is always strongly suppressed. This hints at some common 



underlying physical mechanism which scales with the plasma 
temperature. 

A whole host of explanations have been put forth (for 
reviews, see Fabian et al (2001), Kaastra et al. (2004) 
and references therein). Many center on heating the gas, to 
prevent it from cooling below some minimum temperature 
(e.g. through AGN outflows, electron thermal conduction, 
mergers, or magnetic reconnection) . This suffers from two 
generic problems: firstly, it is notoriously difficult to heat 
gas in a stable manner; for instance, conduction tends to 
either make the gas isothermal (contrary to observations), 
or allow cooling to continue unimpeded, unless the conduc- 
tion suppression factor is excessively fine-tuned (Bregman 
& David 1988; Meiskin 1988) (though see Kim & Narayan 
(2003) for a countervailing claim). Secondly, such solutions 
are much more effective at altering the normalization of the 
emission measure (by reducing M), rather than its shape. 
Indeed, for consistency with observations, heating must ei- 
ther completely suppress local thermal stability or heat the 
gas all the way back up to the ambient temperature. If heat- 
ing delays but does not halt cooling at some lower tempera- 
ture, this will increase the emission measure of low temper- 
ature gas over predictions of the standard isobaric cooling 
model, rather than reduce it. Another possibility is to cool 
gas rapidly by some process which does not involve X-ray 
emission (e.g. turbulent mixing, dust cooling), though such 
proposals to date suffer from similar difficulties in that it 
is difficult to naturally identify a characteristic temperature 
scale at which such effects become important (as is required 
by observations). Other proposals (more similar in spirit to 
the route taken in this paper) seek to suppress the emis- 
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sion lines by resonant scattering or absorption (Peterson et 
al 2001), or by an inhomogeneous metallicity distribution 
which significantly alters the predicted emission line spec- 
trum (Morris & Fabian 2003). However, none have proven 
wholly successful, and the problem remains unsolved. 

In this paper we pursue another avenue: breaking the 
expected correspondence between the thermal state of cool- 
ing gas and its emission spectrum. If the ionization stage of 
an atom does not correspond to its value in collisional ion- 
ization equilibrium, then the emission line spectrum ceases 
to be a good thermometer of the plasma temperature. From 
equation Q, the luminosity in an emission line at frequency 
v from an ionization stage i is given in the standard isobaric 
cooling flow model (where a — 0) by: 



L V = %- 
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Note the inverse weighting by the cooling function: emis- 
sion at lower temperatures is given progressively less weight 
as the gas cools rapidly through that phase. Since dL v oc 
[xi(T)/A(T)] dT (where Xi is the ionization fraction of stage 
i), if we can delay the appearance of an ionization stage i 
until low temperatures when cooling is rapid, the emission 
measure of its associated lines will be severely reduced. This 
could happen, for instance, if the gas falls out of collisional 
equilibrium; however, for conditions in cooling flow clusters 
tree <C icooi, so this is unlikely. On the other hand, the gas 
can also be overionized if there is some additional source 
of ionization. We consider two sources: collisional ionization 
from a non-Maxwellian tail of suprathermal electrons, and 
photo-ionization by free-free emission from hot intracluster 
gas. Both of these effects are associated with the large reser- 
voir of hot gas in the cluster, and the associated ionization 
scales naturally with the temperature of the ambient gas. 
Thus, the onset of excess ionization takes place at some 
fraction of the cluster virial temperature. This would natu- 
rally explain the universality of the low temperature cut-off 
and the reason why it always appears at some fixed fraction 
(~ 1/3) of the temperature of the hot gas, features which 
none of the presently popular models can explain. For in- 
stance, it is difficult to understand why heating mechanisms 
like AGN heating should target low temperature gas at the 
expense of high temperature gas, particularly when cool gas 
occupies a very small fraction of the cluster by volume. 

Despite the above promising characteristics, we find 
that excess ionization fails to solve the soft X-ray cooling 
flow problem, for relatively general and model-independent 
reasons. In clusters, suprathermal electrons are much more 
efficient at heating the gas via Coulomb collisions rather 
than ionizing the metals (unlike in the solar corona, 
where their ionization effects have been invoked to explain 
anomolous ionization patterns, e.g. Scudder (1994)). Pho- 
toionization can certainly produce excess ionization at low 
temperatures with little associated photo-heating, since the 
the ionizing photons only target the metals, which are a rare 
tracer species. However, the same excess ionization reduces 
the efficiency of gas cooling at low temperatures which ul- 
timately increases the emission measure of low ionization 
stage emission lines. This arises because the missing emis- 



sion lines are produced at temperatures where line cooling, 
rather than bremstrahhlung, is the dominant cooling mecha- 
nism; one therefore cannot alter the ionization properties of 
the gas without altering its cooling properties. Despite the 
failure of these two attempts, the generality and tempera- 
ture scaling of such alterations to standard isobaric cooling 
emission spectrum models is very attractive, and we feel, 
under-explored. It is possible that some variant of these ideas 
may ultimately prove successful. 



2 COLLISIONAL IONIZATION: 
SUPRATHERMAL ELECTRONS 

Let us consider a scenario where the electron velocity distri- 
bution is not strictly Maxwellian, but has a non-Maxwellian 
tail of high velocity electrons. In this case, the usual cor- 
respondence between plasma temperature and metal ioniza- 
tion stage may be broken: while the gas can cool down to low 
temperatures, the metals are maintained at a relatively high 
ionization stage by the non-Maxwellian tail of suprathermal 
electrons. Such ionization effects due to a non-Maxwellian 
tail is thought to take place in the solar corona (Roussell- 
Dupre 1980, Scudder & Olbert 1979, Cranmer 1998), where 
the electron mean free path is comparable to the tempera- 
ture scale height. Furthermore, if the typical energy of elec- 
trons in the suprathermal tail scales with the cluster tem- 
perature, this would explain why the cooling cutoff scales 
with the cluster temperature. Below, we build simple toy 
models in which a small suprathermal tail attenuates the 
line flux from lower metal ionization stages by the required 
amounts. Nonetheless, here we show that this promising sce- 
nario does not hold for typical conditions in clusters: any 
non-Maxwellian tail of electrons which ionizes metals above 
their equilibrium ionization stage at a given temperature 
will cause an unacceptable amount of Coulomb heating. 

2.1 Excess Ionization with a Non-Maxwellian Tail 

The presence of a significant population of suprathermal 
electrons in cluster plasma is usually ignored. This is jus- 
tified on the basis that the mean free path of electrons is 
generally much smaller than the temperature scale height. 
A minimal estimate of the temperature scale height Lt is 
the Field length Xf, which can be obtained in a local sta- 
bility analysis by assuming balance between conduction and 
cooling (Field 1965): 



\f = 



f^T 1 ' 2 
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= 150/ 1/2 T 7 3/2 n:J(lnA c )/40r 1/2 kpc(3) 



where lnA c = 29.7 + ln(n~ lll2 Te) is the Coulomb logarithm, 
ks is the Spitzer coefficient for thermal conduction, and / 
is some unknown parameter which accounts for the reduc- 
tion of the electron mean free path (and thus the thermal 
conductivity) due to tangled magnetic fields or plasma in- 
stabilities. By contrast, the mean free path of electrons due 
to Coulomb collisions is: 
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Since A e <C Af, it would appear that the assumption of clas- 
sical conductivity is valid, and further that suprathermal 
electrons cannot play an important role in cooling flows. 
However, a non-Maxwellian tail of suprathermal electrons 
could plausibly arise in two cases: (i) plasma instabilities 
such as whistlers (Pistnner, Levinson & Eichler 1996) scat- 
ter electrons below some energy -Ecrit, but are relatively in- 
effective at scattering electrons about this energy. Thermal 
conduction will be significantly reduced if most of the heat- 
bearing electrons at the peak of the Maxwellian distribution 
have significaantly attenuated mean-free paths. The Field 
length Xf oc 1// and thus the temperature scale height 
Lt will therefore be severely reduced (indeed, sharp tem- 
perature drops indicating conduction suppression have been 
observed in clusters (Ettori & Fabian 2000)). Nonetheless, 
suprathermal electrons with A e > Lt can penetrate the cool- 
ing gas and ionize the metals. Indeed, Chandra has observed 
a sharp temperature drop in gas from 10 to 5 keV on scales 
of order ~ 10 — 15kpc (Markevitch et al. 2000), comparable 
or less than the electron mean free path; such conduction 
suppression processes may be responsible for the survival 
of such sharp fronts, (ii) Jets or shocks in the intraclus- 
ter medium introduce a secondary population of suprather- 
mal/relativistic electrons. Indeed, the observation of syn- 
chrotron radiation from clusters indicates that such a pop- 
ulation of electrons must exist (e.g., Sarazin (1999)). 

To demonstrate the effect of such a suprathermal tail 
on ionization balance, we build a simple toy model. We con- 
sider an isobarically cooling filament of gas embedded in a 
the ambient ICM at the cluster virial temperature of 8 keV, 
with initial density n ~ 10 _3 cm -3 . We assume a plane- 
parallel geometry. Our model has one free parameter, the 
size of the filament H, which we parameterize as a = H/X, 
where A ~ 20kpc is the mean free path of an 8 keV electron 
in the ambient ICM. We solve the Fokker-Planck equation 
for the electron distribution function, taking into account 
the effect of Coulomb collisions. Our approach is similar to 
that of Owocki & Canfield (1986) in the solar corona. We 
do not attempt to solve the full Fokker-Planck equation (e.g. 
Shoub (1983)), a computationally very demanding task. In- 
stead, we use a linearized BGK (Bhatnagar, Gross & Krook 
1954) approach in which the collision operator is given by a 
phenomenological relaxation term. The collisional thermal- 
ization rate is simply given by the rate at which high- velocity 
test electrons are deflected by 90°. While in a strict sense 
the BGK method does not conserve energy or momentum, 
in practice it agrees well with more elaborate calculations 
(Ljepojevic & MacNeice 1988), and both the great ease of 
caculation and the physical insight provided recommends it 
for our use. We ignore the effects of a polarizing electric 
field (e.g., see Bandiera & Chen 1994). With these approxi- 
mations, we solve: 

df 

where / is the electron velocity distribution function we wish 
to determine, f*(v, z) is a Maxwellian distribution with tem- 
perature T(z), v(r,v) = 167re 4 n(r)lnA/(m^u 3 ) is the colli- 
sional thermalization rate, and /i = v z /v. This equation is 
strikingly reminiscent of the equation of radiative transfer 




E (keV) 

Figure 1. The computed electron distribution functions in the 
cooling clump at 1 keV for different a = H/X; also shown are 
Maxwellian distributions for gas at 1 keV an d 8 keV (the tem- 
perature of the ambient ICM). At very high energies, the clump 
is transparent to the hot electrons and the distribution func- 
tion joins smoothly over to that of the ambient ICM. The non- 
Maxwcllian tail is progressively stronger for smaller clumps (i.e. 
smaller a), which have lower optical depth to Coulomb scattering. 

if^- — —I u + S v , which has the formal solution I v (t v ) = 
I„{Q)e~ T " + P"" e" (T "~ T - ) S' I ,(r^)dr,'. Similar, the solution to 
equation © consists of a weighted sum of Maxwellians at 
different temperature, with the weighting given by the at- 
tenuation exp(— r) due to Coulomb shielding. We show the 
computed distribution functions for a — 0.5, 1, 2 and fi = 1 
in Figure At very high energies, the clump is transpar- 
ent to the hot electrons and the distribution function joins 
smoothly over to that of the ambient ICM. 

Given the new distribution function f(v), we compute 
the new ionization and recombination rates for Fe ions in the 
plasma. For a collisional process with cross section a(E), the 
rate coefficient is given by: 

r=( — ) / x 1/2 a(xkT)f(x)dx (6) 
V me 1 J. th 

where x t h = E t h/kT, and E t h is the threshold energy for 
a particular process. We obtain the cross sections for all 26 
ionization stages of Fe for collisional ionization (both direct 
and excitation autoionization) and recombination (both ra- 
diative and dielectronic) from parametric fits by Arnaud & 
Raymond (1992). While the recombination rates are only 
mildly affected by the non-Maxwellian tail of high energy 
electrons, ionization rates are very significantly affected, and 
can be orders of magnitude larger: comparable to the ion- 
ization rate in the ambient medium, for a cooling clump suf- 
ficiently transparent to suprathermal electrons. If we then 
assume ionization equilibrium-a good assumption since the 
recombination time is very short compared to the cooling 
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Figure 2. Bottom "panel: The fractional abundance of Fe XXIV 
(whose emission lines are observed, dotted lines) and Fe XVII 
(which emission lines are not observed, dashed lines) as a function 
of clump temperature. Shown as solid lines are the abundances for 
purely Maxwellian electron DFs. For small clumps a < 0.5, the 
abundance of Fe XVII is strongly suppressed until very low tem- 
peratures, implying that Fe XVII emission lines will not be seen. 
Top panel: The mean charge of Fe as a function of clump temper- 
ature, for clumps of different sizes. Shown in black is the expected 
evolution for a purely Maxwellian electron DF. For small clumps, 
suprathermal electrons keep Fe ions at high ionization stages even 
at low temperatures, reducing the usefulness of observed lines as 
a temperature diagonostic. 



time, e.g. for Fe XVII, t rec ~ 10 6 nl2yr < t coo i ~ lC^nl^yr- 
we can then compute the local ionic abundance Xi for all 
ionization stages as a function of clump temperature. In 
Figure [5] we display the fractional abundance as a function 
of local temperature for selected ionization stages, as well 
as the mean ionization stage as a function of temperature. 
Compared to a Maxwellian DF, the mean ionization stage 
is significantly higher, particularly at low temperatures. In 
particular, the abundance of lower ionization stages as Fe 
XVII is severely suppressed until low temperatures, when 
the cooling time is very short. 

We can quantify the suppression of low ionization stage 
cooling lines. Standard isobaric cooling flow prescriptions 
have an emission spectrum given by equation indepen- 
dent of geometry. The inverse weighting by the cooling func- 
tion implies that little emission is associated with gas cooling 
at low temperatures, when the cooling time is short. This 
allows us to compute the suppression factor of line emission 
from a particular ionization stage i: 



fT 
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A(T) 
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0.001 




0.01 - 
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Figure 3. Emission line suppression factor y su PP ress f or eacn Fe 
ionization state, as given by equation {7J, for a = 0.5,1,2. A 
purely Maxwellian electron distribution function corresponds to 
/■suppress = 1. For distribution functions with a non-Maxwellian 
tail, line emission from the higher ionization stages is enhanced, 
while emission from the lower ionization stages is strongly sup- 
pressed. For a = 0.5, 1, the suppression far exceeds the minimal 
amount required to explain observations (j= u PP rcss < q 3 f or p e 
XVII and lower in this case). 



where 2;^ upra an( j x ^ 



are the fractional abundances 



A(T) 



of ions assuming a distribution function with a tail of 
suprathermal electrons, and a locally Maxwellian DF respec- 
tively. We have assumed that the local cooling rate A(T) 
is not significantly affected by suprathermal electrons (we 
shall see in il4.3l that this assumption is not strictly correct, 
but here the suppression of the lower ionization stages is 
so strong that this correction is unimportant). We there- 
fore adopt a fit to the cooling function of Sutherland & 
Dopita (1993), for Z = O.3Z . We show the results in Fig- 
ure lO, which show emission from lower ionization stages 
to be severely suppressed. By comparison, for example, the 
XMM observations of Abell 1795 which show a lack of Fe 
XVII emission lines require that their emission measure be 
suppressed by at least a factor of 3 (Tamura et al. 2001). 
Hence, in principle such effects could easily explain the lack 
of emission lines associated with low temperature gas. 

Note that the suprathermal tail of the electron dis- 
tribution function in cooling gas is simply given by the 
high-energy Maxwellian tail of the ambient cluster medium. 
There is therefore be a correlation between the virial temper- 
ature of the cluster and ionization stage at which emission 
becomes strongly suppressed. The cutoff temperature below 
which line emission is suppressed is always some fraction of 
the cluster virial temperature. This self-similar feature of the 
excess ionization model, which closely mirrors what is seen 
in observations, is a very promising characteristic. It does 
not emerge naturally in most standard heating models. 
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2.2 Unacceptable Coulomb Heating 

Thus far, we have not self-consistently included the heating 
due to suprathermal electrons on the temperature profile 
of the clump. For simplicity let us consider a suprathermal 
population of electrons of temperature Th and density rih 
penetrating a cooling clump of gas at temperature T c and 
density n c . Note that Th and rih do not necessarily corre- 
spond to the temperature and density of the ambient cluster 
gas; they are determined by the details of the electron scat- 
tering process. In general T h > Tdustor and n h < n c i ust er, 
for high energy electrons far in the Maxwellian tail. We al- 
low ourselves arbitrary freedom in the two variables Uh,Th 
to see if there exists any form of the suprathermal tail which 
would ionize the metals, while still allowing the gas to cool. 

The rate of energy loss of a highly suprathermal electron 
Th 2> T c is independent of T c and given by (Schunk & Hays 
1971): 

dT 4 7 rn c e 4 1 / m s vl \ 
at m e v e Kje^LOp J 

where uj p = 4nn c e 2 /m e is the plasma frequency. This yields 
a net Coulomb heating rate: 

n c T » 5.5 x 10" 14 nhn c r h " 1/2 ergs _1 cm" 3 (9) 

The ionization timescale for the ionization stage i is given by 
t ion = n l /hl on = l/n h Oi{v h )v h (where cn{v h ) is the velocity- 
dependent ionization cross-section of ionization stage i), 
while the recombination time is given by t rec — l/a rec n c , 
where a rec — c^+i + af+i is the sum of the radiative and di- 
electronic recombination rates. The ionization cross-section 
and recombination rates for Fe are evaluated from paramet- 
ric fits by Arnaud & Raymond (1992). We also define a heat- 
ing timescale thcat = kBn c T c /n c T and a cooling timescale 
tcooi = k B n c T c /n 2 c A(T c ). 

For definiteness, we consider the effect of suprathermal 
electrons on Fe XVII, which is not seen in XMM observations 
of cooling flow clusters (Peterson et al 2001). We need only 
consider the temperature range T c w 0.2—0.8 keV, where the 
abundance of Fe XVII becomes appreciable in a Maxwellian 
plasma. Since t rcc ,t coo \ oc l/n c , we find: 

i rcc » 0.3 - 1 X 10~ 3 i coo i; 0.2keV <T C < 0.8keV (10) 

independent of n c . Likewise, since thcat, tion oc 1/uh, we find: 

tion - 10t hcat \ (11) 

independent of rih,Th. The independence from Th arises be- 
cause tionAheat oc a(v h )T h , and at high energies the Bethe 
approximation a oc v^ 2 holds. For the abundance of FeXVII 
to be strongly suppressed, we require U on < t Icc . However, 
from equations (1101 and I1H . this implies 

thcat < 10~ 5 - W~ 4 t CO oi (12) 

for the relevant temperature range T c « 0.2 — 0.8 keV, 
and independent of n c ,n,h,Th. Thus, any population of 
suprathermal electrons which is sufficient to ionize Fe 
XVII at low temperatures would cause unacceptably large 
amounts of Coulomb heating which cannot be radiated away, 
and the cooling clump would be rapidly evaporated. Such a 



large suprathermal tail would not permit gas to cool signifi- 
cantly below the ambient cluster temperature to begin with. 
This is a robust conclusion despite our order-of-magnitude 
approach: the heating rate implied by the required ionization 
exceeds the gas cooling rate by 4-5 orders of magnitude. 

Our result can be understood from the fact that the 
cross-section for inelastic collisions of electrons with atoms 
is much smaller than the cross-section for elastic Coulomb 
collisions with other electrons by at least the value of the fine 
structure constant, ~ 1/137. The effectiveness of frequent 
small-angle encounters at small impact parameters decreases 
the electron-electron equilibration timescale by another 2 or- 
ders of magnitude (Spitzer 1978). Thus, the plasma is heated 
much more rapidly by suprathermal electrons than it is ion- 
ized. 1 

In summary, the primary effect of a non-Maxwellian tail 
in the electron velocity distribution in cluster gas is on the 
plasma temperature rather than the ionization stage (ex- 
cept possibly for very low ionization stages, which have long 
recombination times). Thus, the ionization stage of Fe re- 
mains a good diagnostic of plasma temperature, and the 
lack of emission lines from low ionization stages seems to 
suggest that the gas cannot cool below some cutoff temper- 
ature T cr it ■ 



3 EFFECTS OF NON-THERMAL PRESSURE 
SUPPORT 

In this section, we very briefly consider the effects of a source 
of non-thermal pressure support in cooling gas, either in 
the form of magnetic field or cosmic rays. This will be rel- 
evant for evaluation of gas densities in our consideration of 
photoionization in ^\ Observations of Faraday rotation in 
the cores of clusters indicate the presence of highly chaotic 
> 5 — 10/iG magnetic fields (Clarke, Kronberg & Bohringer 

2000) , which in some cases reach equipartition values (En- 
blin et al 1997). In addition, the observation of synchrotron 
emission, jet and radio galaxy activity mean that the energy 
density in magnetic fields and cosmic rays could be consid- 
erable, and estimates show that it can be comparable to 
the thermal energy density of the gas (Enblin et al 1997; 
Volk, Aharonian & Breitswerdt 1996). Finally, cosmic rays 
could also be accelerated in gravitational shocks due to ei- 
ther accretion (Miniati et al 2001) or mergers (Donnelly et al 

2001) . As the gas cools, the non-thermal component evolves 
adiabatically and can eventually dominate the overall pres- 
sure support. Since the gas no longer cools isobarically, this 
has two important consequences for us: (i) the emission line 
spectrum differs from that of standard isobaric cooling mod- 
els; (ii) gas at low temperatures cools at lower densities, re- 
ducing the recombination rate and allowing it to be more 
easily photoionized (see [Q). In the limit where non-thermal 

1 For similar reasons, Canizares et al (1993) concluded that for 
a given mass flux, the amount of X-ray emission from a 'heating 
flow' (where cool gas is evaporated and implusivcly heated to the 
ambient cluster temperature) is ~ 10 — 4 times the emission from 
a cooling flow. 
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pressure dominates, cooling becomes isochoric rather than 
isobaric, leading to a factor of at most ~ (5/2)/(3/2) ~ 5/3 
reduction in the strength of low temperature emission line 
features; thus effect (i) by itself cannot account for the fac- 
tor > 3 reduction in the strength of low temperature emis- 
sion lines. However, effect (ii) implies that the recombination 
times could be longer by up to a factor ~ (Ti/T) from the 
isobaric cooling case (where Ti is the initial background tem- 
perature, and T is the temperature at which a particular 
ionic stage becomes abundant). This permits photoioniza- 
tion to dominate over recombination, and the latter effect 
could be quite important. Here we only consider the dynam- 
ical effects of a non-thermal component; we shall consider 
the possible implications of cosmic-ray heating or magnetic 
reconnection (when the non-thermal component does not 
evolve adiabatically) in a subsequent paper. 

What is the effect of non-thermal pressure support on 
a thermally unstable clump of gas? If we model the non- 
thermal pressure as a relativistic gas with adiabatic index 7 
(henceforth we shall assume 7 = 4/3) then the total pres- 
sure: 



nksT + Knttl 1 = const 



(13) 



is conserved during the cooling process, if the cooling time 
is much shorter than the characteristic flow time, t coo i <C 
tfiow ~ R/u, which is generally the case (we have also as- 
sumed that the magnetic field lines and any associated cos- 
mic rays are frozen into the plasma). Given the initial gas 
pressure, non-thermal pressure and temperature of the gas, 
P g as,i, PNT,i,Ti, we can solve equation 1131 for the density as 
a funtion of temperature n(T), assuming that the adiabatic 
constant Knt = PcR,i/nJ is conserved. The gas density at 
given temperature is always less than that assumed by the 
isobaric cooling model n(T) < rii SO baric = riiTi/T; also, at 
low temperatures non-thermal pressure support eventually 
dominates, and the gas cools almost isochorically. These ef- 
fects are illustrated in Figure 21 for different values of the 
initial non-thermal pressure /; = PNT,i/Ptot,i- For cluster 
cores, it is reasonable to assume up to equipartition values 
(fi ~ 0.5) for the non-thermal pressure, before the onset of 
thermal instability. We will use equation 1131 to compute 
the densities n(T) when calculating recombination rates in 
□ 

We can easily compute the effect of the non-thermal 
pressure support on the emission measure of the gas as a 
function of temperature. It now becomes: 



dEM 
dT 



(T 2 ) = 



3 712/ni 
2 + T1/T2 



1 



M kr 



1 J p,m v A(T) ' 



(14) 



where m , Ti are the initial temperature and density, and 
n^iT^) is obtained from equation 1131 In the limit where 
non-thermal pressure is unimportant, 712/ni = T1/T2, and 
the prefactor — > 5/2, for isobaric cooling; when non-thermal 
pressure completely dominates, n-z w m and the prefactor 
— > 3/2, for isochoric cooling. 

Non-thermal pressure support will also retard gas in- 
flow and reduce the amount of work done on the cooling gas 
by the gravitational potential. Note that the compression 
of gas in the dark matter potential will make little differ- 



o 



o 



>. 




-isobaric fi— 0.1 




- isochoric 



-em <r 



em «r 



0.2 0.4 0.6 0.8 1 

T/T 

Figure 4. Bottom panel: The emission measure as a function of 
temperature for models with initial fractional non-thermal pres- 
sure fi = PNr/Ptot at temperature T . Nonthermal pressure 
support results in smooth interpolation between the isobaric and 
isochoric cooling regimes and cannot account for the observed 
emission measure, where EM oc T a , with a sa 1 — 2. Top panel: 
Density as a function of temperature, normalized to the density 
the gas would have if it cooled isobarically. Non-thermal pres- 
sure become progressively more important at lower temperature, 
leading to significantly lower densities than if the gas cools iso- 
barically. This could increase recombination times and delay the 
formation of low ionization stage ions, if there is a significant 
photoionizing flux (see Q . 



ence to the overall gas cooling rate, since for temperatures 
and densities characteristic of cooling flow clusters the gas 
cooling time depends almost exclusively on its entropy. In 
particular, at fixed entropy the cooling time is almost inde- 
pendent of the gas density; thus, the total cooled gas mass 
is almost unchanged whether the gas is adiabatically com- 
pressed by the gravitational potential or not (Oh & Benson 
2003). Thus, the presence or absence of non-thermal pres- 
sure support makes little difference in terms of the actual 
mass drop-out rate (although it can have a large effect on 
the inferred mass drop-out rate, since the bolometric X-ray 
luminosity of the cluster can be significantly reduced when 
non-thermal pressure restricts inflow). However, its main ef- 
fect is in altering the normalization rather than the shape 
of the X-ray spectrum (Peterson et al 2002; Nulsen 1998; 
Allen 2000). In particular, compression by the gravitational 
potential makes little difference to lines which appear at low 
temperatures, since t coo i -C tfi ow » R/c s , and the gas cools 
before it travels a significant distance. As noted by Peter- 
son et al (2002), even if there were no mass drop-out, the 
maximal effect of gravitational inflow would be to replace 
in equation (when a — 0) the factor 5/2 — * 3/2 — \t/\ p , 
where At is the radial logarithmic temperature gradient and 
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X p is the radial logarithmic density gradient; At/A p is always 
observed to be negative in the cores of clusters. 



4 PHOTOIONIZATION 

4.1 Is Photoionization of Metals in Clusters 
Plausible? 

From the preceding discussion in H2.2I it is clear that if we 
wish to break the correspondence between metal ionization 
stage and plasma temperature by ionizing the metals, we 
need to use an ionizing source which injects little heat per 
ionization. Photoionization is ideal for this purpose: it only 
targets the metal ions (the Compton optical depth r < 10~ 2 
through clusters is small); and since metal ions are a rare 
tracer species, photoheating is generally unimportant. Fur- 
thermore, in clusters a large reservoir of hot X-ray emitting 
gas is available to provide ionizing photons. In clusters it is 
generally assumed that the plasma is in the coronal limit 
and that the ambient radiation field is too weak to play a 
role in the ionization properties of the gas. Below we per- 
form a simple order of magnitude estimate to show that in 
fact photoionization effects may be non-negligible; indeed, it 
has been proposed (Donahue & Voit 1991; Voit, Donahue & 
Slavin 1994) that the optical filaments seen in clusters are 
powered by self-irradiation by the surrounding hot cluster 
gas. Once again, since the free-free luminosity scales with 
the virial temperature of the cluster, this would provide an 
explanation of the universal scaling of the cut-off temper- 
ature with the cluster temperature: the harder radiation 
field in hotter clusters can induce photoionization of pro- 
gressively higher ionization stages. For definiteness, we once 
again consider the ionization of Fe XVII, an ion indicative 
of plasma cooling down to 0.2-0.8 keV, which is notoriously 
absent in the spectra of cooling flow clusters (Peterson et 
al 2001). In steady state, the energy advected into a cool- 
ing region should balance the radiative energy emerging. As 
Voit, Donahue & Slavin (1994) argue, the former should be 
~ u(r)v(r), where u(r) is the effective energy density at r 
and v(r) is the inflow velocity. If the cooling occurs isobar- 
ically, then u(r) ps 5P(r)/2. The inflow speed should be of 
order the sound speed in the gas, v(r) ~ c s ~ 300kms~ , 
where the latter is the isothermal sound speed of gas at 
~ IkeV. This energy flux ~ u(r)v(r) is a reasonable esti- 
mate of the ambient radiation field in a cluster. If it mostly 
emerges in the form of free-free emission, the rate at which 
a metal atom is photoionized is: 

Rion « \-j^f iCJi ( 15 ) 

where Ei is the energy of the ionization edge, Oi is the pho- 
toionization cross-section at ~ Ei, and fa is a correction 
factor for the fraction of the flux emerging above Ei\ for a 
free- free spectrum typically ft is of order unity. By contrast, 
the recombination rate is: 

R Iec = a(T)n e = a(T) (JL^ (16) 

where a(T) is the temperature dependent recombination co- 
efficient, and we have assumed that cooling occurs isobari- 



cally. The ratio of the photoionization rate to recombination 
rate is therefore: 

Photoionization ^ k B T ( c B a r \ /k B T\ ( a(k B T)\ 

Recombination ET \a(T) J ~ ' V Ei J \ a(Ei) j 

where the quantities a ~ 10 _11 cm 3 s _1 , a,: ~ 10 -19 cm -2 
have been evaluated for k B T ~ Ei ~ 1.26keV, as is ap- 
propriate for Fe XVII. It is interesting to note that this 
result is independent of the pressure and hence the density 
of the gas: for isobaric conditions, higher densities result in 
a higher recombination rate but also a higher photoioniza- 
tion rate, since the surrounding gas has a higher emissivity. 
The close match of the photoionization and recombination 
rates in this rough estimate warrants further investigation. 
In addition, we shall see that other sources of photoionizing 
radiation such as AGN or X-ray binaries could be even more 
important. 

4.2 Estimating the Ambient Radiation Field 

Let us quantify the level of photoionization needed to 
explain the observational results. We define a critical 
flux F cr it, such that the photoionization rate Ri 0n = 
J dEF cr u(E)ai(E) is roughly equal to the recombination 
rate to Fe XVII, R r ec(T c ) = a(T c )n e , c where n c is the elec- 
tron number density at T c = 0.4keV, the temperature at 
which the abundance of the Fe XVII ion peaks. For free- 
free emission from a plasma at temperature Th ~ 5keV, the 
critical flux at the ionization edge Ei = 1.26keV of Fe XVII 
is: 

FcritiEi) = 4.4 x 10 s ( n ™ c \ keVcm-V'keV- 1 (18) 
VO.lcm - ' 5 / 

Due to the flat spectrum of free-free emission up to ~ 
Th, this depends only very weakly on the temperature of 
the hot plasma Th, as long as Th ^ Ei = 1.4keV. For 
a power law ionizing spectrum F oc E~ r where V — 
0.7, we obtain a very similar result, F^f t N (Ei) — 4.2 x 
lO^nc/O.lcm-^keVcm-V^eV" 1 . This gives us a rough 
estimate of when photoionization is non-negligible; hence- 
forth, we will normalize all our results to F cr it{Ei). 

Let us first consider the effects of free-free emission from 
hot cluster gas. We estimate the internal radiation field of 
a cluster following Sazonov, Sunyaev & Cramphorn (2002) , 
who did so in order to constrain resonant X-ray scattering. 
Consider a point some distance r from the center of a cluster. 
Integrating over all solid angle, the total flux of radiation 
received at that point is: 

dfj, J dr'e{R) (19) 

where e(r) is the emissivity profile of the gas, and R 2 = 
r 2 + r' 2 — 2rr' fi. The continuum emissivity of cluster gas in 
bremsstrahlung is: 

e(E) = 2.3 x W' 20 T' 1/2 exp{~E/k B T)n 2 g B (T,E) (20) 

erg cm s kev 

where g B (T,E) w (E/k B T)~ 0A is the Gaunt factor. Let us 
consider a cluster with a beta-law radial density profile: 
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(l + r 2 /r§)3/3/2 



(21) 



where n is the central electron density and r c is the core ra- 
dius. Equation 1191 can be integrated numerically to obtain 
the X-ray flux at any point in the cluster. However, there are 
two interesting analytic limits. At the center of the cluster: 



F(r = Q,E) = 10 



7 r( 3/3-l/2) 
1X3/3) 



/ n 
\ 0.1 cm 



/ k B T \ 
UkeVy 



-1/2 



x °' 4 exp(— x) keVcm 3 s 



E 



\ lOOkpc 



VO.lcm 



lOOkpc 
^eV" 1 

n c /n 



— j exp[-x + Xi] 

where x = EjksT . Also, for the special case of j3 = 2/3, the 
flux traces the density profile exactly: 



F(r,E)=F(0,E) 



1 + (r/r c ) 2 



oc n(r) 



(23) 



Since F cr i t oc n, this implies that F/F CI it = const through- 
out the cluster, i.e., the decreasing flux and photoionization 
rate at the outer parts of the cluster gas is tracked exactly 
by the decreasing densities and hence decreasing recombi- 
nation rate. For gas in an isothermal potential, this implies 
that photoionization suppression of lower ionization stages 
is constant throughout the cluster. Note that, even if the as- 
sumption of isothermality is dropped, F oc n is a very good 
approximation (since the emissivity depends largely on den- 
sity and is only weakly temperature dependent). This could 
provide a natural explanation as to why the cluster gas ap- 
pears locally isothermal at each radius (Molendi & Pizzo- 
lato 2001; Matshushita et al 2002): even if the cluster gas 
is multi-phase, the cooler gas is overionized for its electron 
temperature, and thus the emission measure of lower ioniza- 
tion stage ions indicative of cool gas is unmeasurably small. 
Furthermore, this effect is relatively independent of cluster 
radius. 

Other possible sources of X-ray ionizing photons are 
AGN and X-ray binaries. The flux from an AGN is Fagn = 
T E 1 ^ a , or numerically: 



47t(E. 



(2-a)L x 
T~ 



Fagn = 1.2 x 10 s 



/ M BH \ ( L x /L 



(—) 

VlkeV/ 



^ 1O 9 M 

-0.7 



0.05 



V lOkpc 



-2 

24) 



keVcm 3 s J keV 



where a ~ 1.7 is the photon index, Lx is the AGN lumi- 
nosity in the energy range E\ — E2 ~ 1 — lOkeV range, and 
we have assumed that ~ 5% of the bolometric luminosity 
emerges in this range (Elvis et al 1994). Assuming again a 
beta-law profile for the density, we find that: 



Fagn(E) « 50F clit (Ei 



on 



\ lOOkpc 



+ 



1O 9 M 
2" 



lOkpc 



Lx/Ly 
O05 

n c /n\ ( E 



/ n 
VO.lcm 



VlkeV J 



It is therefore apparent that particularly near the center 
(where soft X-ray emission lines are most egregiously miss- 
ing), an X-ray bright AGN could produce photo- ionizing 
radiation orders of magnitude higher than required to dom- 
inate collisional processes. It is also interesting to note that 
since Fagn oc r~ 2 and F cr it oc n oc r~ 2 , outside the cen- 
tral density core the AGN flux is also always some con- 
stant multiple of the critical flux, similar to the free-free 
radiation case. The invocation of AGN to provide ioniz- 
ing radiation rather than the ambient hot cluster gas is 
somewhat less desirable, because such sources have a short 
duty cycle and are not seen in all cooling flow clusters (al- 
though Thompson and resonant line scattering could trap 
their radiation for some time and boost the X-ray surface 
brightness by a factor ~ 3 — 10 over the free-free continuum 
(Sazonov, Sunyaev & Cramphorn 2002). Furthermore, the 
universal scaling of the cut-off temperature with the ambient 
plasma temperature is lost, although one could argue that 



Fagn oc M bh oc a CD oc T 



,/2 



where cjcd is the velocity 



dispersion of the CD galaxy and a is the slope of the Mbh—<? 
relation. Nonetheless, it is clear that photoionization effects 
could be important. In the following section we consider if a 
photo-ionizing flux of arbitrary magnitude could be respon- 
sible for the observed emission line spectrum. 



4.3 Ionic Abundances and Line Suppression 

Let us now perform a more careful calculation of the effects 
of photoionization on the emission spectrum of cooling gas. 
Since t ICC <C tcool, the assumption that the metal ionization 
stages assume their equilibrium values (where photoioniza- 
tion and collisional ionization balance recombinations) is a 
good one. However, when computing the emission line spec- 
tra, several effects need to be taken into account. One is that 
the photoionization cross section take into account Auger 
multi-electron ejection. This process increases the impact of 
photoionization of inner shell electrons (as many as 8 elec- 
trons can be removed from the photoionization of the Is 
shell), and couples non-adjacent stages of ionization, which 
makes it necessary to iterate on the ionization solution. The 
second is the change in the cooling function due to the de- 
viation of the ionic abundances from their collisional equi- 
librium values. 

We therefore run CLOUDY (Ferland et al 1998), which 
can compute both of these effects. We present the results of 
these calculations in Figs and |S| In Figure we show 
the result of photoionization on isobarically cooling gas, 
for different values of the radiation field. Only large values 
~ 100-F cr it of the radiation field have a significant impact. 
Furthermore, the net effect of the X-ray background is sim- 
ply to increase the mean charge of the system at a given 
temperature, and shift the appearance of lower ionization 
stages to lower temperatures. Naively, one might have ex- 
pected the X-ray background to become increasingly im- 
_Wrj;ant at low temperatures, since lower ionization stages 
3 nave lower ionization thresholds (and thus are more eas- 
ily ionized). However, the greater densities and recombina- 
tion rates at lower temperatures in isobarically cooling gas 
roughly offset the increase in the photoionization rate. Note 
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Figure 5. The ionic fraction, mean charge, cooling function, and 
cooling-function weighted cmissivity for the FeXVII 15A line and 
FeXXIV llA line for isobarically cooling gas subject to a photo- 
ionizing radiation field Faqn oc u~ 0,7 of strength (dark solid 
line), .F cr it (indistinguishable from no radiation), 10F cr i t (light 
solid line) and 100F cr ; t (dot-dashed line). Only at very high inten- 
sities does the photoionizing background have significant impact. 
The overall effect of an ionizing background is to increase the 
mean charge of a system at all temperatures. Although there are 
more soft photons able to photoionize lower ionization stages, at 
lower temperatures, densities and hence recombination rates are 
much higher, and the two effects roughly cancel. The results of 
calculations with a free-free ionizing background with equivalent 
values of F cr ; t do not differ significantly. 



that a photoionizing radiation field decreases the amplitude 
of the cooling function at low temperatures. This can be 
easily understood: because the metal ions are over-ionized 
with respect to their purely collisional values, the relatively 
cold electrons cannot excite atomic transitions which result 
in cooling radiation. In the limit where the photoioniza- 
tion is very large and this effect dominates, then line cool- 
ing is no longer important and cooling is entirely due to 
bremmstrahhlung, A(T) oc T 1 ^ 2 . We have checked directly 
that photo-heating is entirely negligible in this temperature 
range, even for large (~ 100F cr it) radiation fields. The most 
important quantity is eu ne (T) / A(T) (lower right panel), the 
cooling function weighted emissivity of a given line; the in- 
tegral of this function over all temperatures is directly pro- 
portional to the emission measure of the line, equation 
It is this quantity we wish to suppress in lines associated 
with lower ionization stages such as Fe XVII. Overall, we 
see that for isobarically cooling gas, this quantity is not sig- 
nificantly suppressed, even for high levels of the radiation 
field. Results for a radiation field normalized to the same 
value of F cr it but for a free-free spectrum yield very similar 
results. 

In Figure ©, we show the effects of non-thermal pres- 



Figure 6. The ionic fraction, mean charge, cooling function, and 
cooling-function weighted emissivity for the FeXVII 15A line and 
FeXXIV llA line for gas where only collisional process operate 
and is isobarically cooling (dark solid line) and gas subject to 
a AGN radiation field of strength 10F cr i t which has initial non- 
thermal pressure support Pjvt ~ O.lPtot (light solid line) and 
Pnt ~ 0.5Ptot (dot-dashed line). Because a non-thermal com- 
ponent causes the gas to cool isochorically at low temperature, 
densities and recombination rates are significantly lower at lower 
temperatures, and X-ray photoionization can become significantly 
more important than before. Naively, we would expect from the 
suppression of FeXVII abundance down to lower temperatures 
(upper left panel) that the emission measure of Fe XVII is re- 
duced. However, because of the accompanying reduction in the 
cooling function at low temperature (upper right), the emission 
measure of Fe XVII integrated over all temperatures, is in fact 
increased (lower right). 

sure support, for a AGN ionizing background of 10-F C rit, 
and different levels of initial nonthermal pressure support 
fi = PNr/Ptot- Since the effects of photoionization scale 
as F v /n, if non-thermal pressure support impedes compres- 
sion of the gas to high densities, then (particularly at low 
temperatures) the effects of photoionization could become 
considerably more important. We see that this is indeed the 
case; the appearance of lower ionization stages is suppressed 
at low temperatures. If there were no change in the cooling 
function, then the emission measure of lines associated with 
low ionization stages such as FeXVII would be reduced by 
a factor of a few (depending on the strength of the ioniz- 
ing radiation field), as required to explain the observations. 
However, the reduction of the cooling function at low tem- 
peratures noted above offsets this effect . The net result is 
that the emission measure of lines associated with low ion- 
ization stages is comparable or even larger than the case 
where there is no ionizing radiation. Overall, it is not possi- 
ble to decouple the cooling and ionization properties of the 
gas because the emission lines we wish to suppress appear 
at temperatures when metal line cooling dominates. Thus, 
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this mechanism would only be important in gas where the 
cooling function is independent of the ionization state of 
the metals; for instance, in very low metallicity gas where 
free- free emission cooling dominates. 



5 CONCLUSIONS 

In this paper we have examined two means of decoupling the 
ionization stage of metal ions from the plasma temperature: 
excess ionization from a non-Maxwellian tail of suprather- 
mal electrons, and photo-ionization either from free-free 
emission by hot cluster gas or from internal AGN or X- 
ray binaries. These are very attractive because they should 
be universal and also become important at some fixed frac- 
tion of the ambient plasma temperature, observed properties 
which none of the present-day models based on gas heating 
can explain. Neither of these are successful: suprathermal 
electrons are much more efficient at heating the gas than ion- 
izing the metals, and photo-ionizing radiation decreases gas 
cooling at low temperatures, ultimately increasing the emis- 
sion measure of emission lines associated with low temper- 
ature gas. Nonetheless, they have attractive features which 
may prove useful in elucidating a final solution. In a com- 
panion paper, we shall examine the feasibility of solutions 
involving heating. In particular, we shall examine the global 
stability of cooling flow clusters, addressing the classic cool- 
ing flow problem of mass dropout. 
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